Ferromagnetism in itinerant two-dimensional t 2g systems 
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Motivated by the recent indications of ferromagnetism in transition metal oxide heterostructures, 
we propose a possible mechanism to generate ferromagnetism for itinerant ta s systems in two spatial 
dimensions that does not rely on the coupling between local moments and conduction electrons. We 
particularly emphasize the orbital nature of different bands and show that, when the Fermi level lies 
near the bottom of the upper bands, the quasi-one-dimensional nature of the upper bands leads to 
a large magnetic susceptibility, and thus the ferromagnetic inter-orbital Hund's coupling can easily 
polarize the upper bands and drive the system into a ferromagnetic state. We discuss the connection 
between our mechanism with several itinerant tz g systems that may have ferromagnetic instabilities. 

PACS numbers: 71.30.+h, 71. 10. Ay, 71.45.Lr, 75.30.Fv 



in 

(N 



i 

-f— > 

CO 



I 

C 

O 

o 



> 
00 

in 

o 



Possible ferromagnetism at polar interfaces between 
SrTi0 3 (STO) and other oxides such as LaA10 3 (LAO) 
or GdTi0 3 (GTO)[lHi 

has raised considerable excite- 
ment. Such ferromagnetism is remarkable as the elec- 
trons are believed to reside nearly completely in the t2 g 
bands of the STO, where there are no localized partially 
filled shells to form local moments. Ferromagnetism in 
purely itinerant systems, while envisioned long agoQ, 
is quite rare in practice; most examples may be at least 
partially attributed to local moments formed by partially 
filled d-shells (and often other delocalized electrons). 

In this paper, we discuss the possibility of ferromag- 
netism in t2 g systems of this type. Since ferromagnetic 
polarization inevitably costs kinetic energy, due to the 
Pauli exclusion principle, the key issue is how the kinetic 
energy can be minimized. This may occur when elec- 
trons are close to Mott localization, which requires high 
electron density in the most strict sense, i.e. approaching 
a substantial fraction of unity per atom. Alternatively, 
kinetic energy can be suppressed by Anderson localiza- 
tion, i.e. disorder; for now we concentrate on clean sys- 
tems, however. Otherwise, details of band structure may 
provide reduced dimensionality which also lowers kinetic 
energy. Such a mechanism of inducing ferromagnetism is 
much rarer, but is known to occur. An extreme example 
is the ferromagnetic state at v = 1 in the quantum Hall 
regime, in which a completely flat band is produced by 
the orbital effect of the magnetic field @. 

We would like to make a distinction between ferro- 
magnetism, and correlated electronic phases more gen- 
erally, caused by ioca] or long-range Coulomb interac- 
tions. Most correlated electron physics occurs in transi- 
tion metal compounds and is attributed to the former, 
local Hubbard-like interaction. Since such local inter- 
actions are effective only if there would be a substantial 
probability for uncorrelated electrons to occupy the same 
site, they lead to correlation effects only at high elec- 
tron density. Local interactions of this type are almost 
completely negligible in typical doped semiconductors. 



By contrast, long-range Coulomb interactions become in- 
creasingly important at low carrier density, and dominate 
the interaction physics in semiconductors. When this 
is the case, an effective mass description generally ap- 
plies, and dimensional analysis implies that the impor- 
tant parameter describing the strength of correlations is 
the dimcnsionlcss quantity, r s , which gives the ratio of 
the 1/r Coulomb energy to kinetic energy per electron. 
It increases with decreasing electron density, r s ~ n~ l / d , 
in d dimensions. This "jellium" problem has been ex- 
tensively analyzed by highly accurate quantum Monte 
Carlo simulations, and it is known that quite large val- 
ues of r s are needed to destabilize the ordinary, para- 
magnetic metallic state (r a > 37 ± 5 in two dimensions 
and r s > 75 ± 5 in three dimensions^, [l(|). Since the 
corresponding electron density is extremely small, this 
itinerant ferromagnetism in low density systems due to 
this Stoner's mechanism is extremely fragile if it exists 
at all, and we disregard it in the following. 

We are led therefore to consider high electron den- 
sity ferromagnetism. In the aforementioned LAO/STO 
and GTO/STO interfaces, there is an intrinsic mecha- 
nism for high carrier density: the polar discontinuity (ll| . 
LAO and GTO have a structure of polar (001) layers: 
La 3+ 2_ has a net charge of +1 per unit cell, while 
Al 3+ (0 2 ~)2 has a net charge of —1 per unit cell (the 
same counting holds for GTO). STO by contrast is non- 
polar. At an ideal (i.e. without atomic reconstruction or 
compensating defects) interface between two such mate- 
rials, an electron gas is predicted to arise with a carrier 
density of half an electron per two-dimensional unit cell. 
This translates, using the unit cell of STO, into a two- 
dimensional carrier density of n = 3.5 x 10 14 cm -2 , which 
is extremely large by semiconductor standards. But is it 
high enough for ferromagnetism and correlation physics? 

It is important to realize that a large two-dimensional 
density does not necessarily imply a large electron density 
per atom. Due to the hnite width of the confined wave- 
function of electrons in STO, which has been estimated 
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now in many calculations j!2Hl5j . the electron occupation 
per site should be less than 20%, even for the Ti sites 
nearest the interface. This fact alone makes strong ferro- 
magnetism hard to understand. Modern computational 
studies have put strong restrictions on ferromagnetism in 
Hubbard models 0. The most recent studies argue that 
ferromagnetism is absent in the two-dimensional Hub- 
bard model for fractional site occupation x < 0.7, even 



when the on-site Hubbard interaction U — > oo 

The remainder of the bulk of this paper explores the 
possibility of weak (i.e. small moment, low temperature) 
ferromagnetism in t2 ff systems at carrier densities which 
are high by semiconductor standards, but still achievable 
near STO interfaces and in bulk doped STO. In particu- 
lar, we will assume that the fractional filling of all bands 
is sufficiently small that only states near the bottom of 
a band need to be considered, but that the contribution 
from on-site interactions is not negligible. We will refer to 
this as intermediate-high carrier density. We first discuss 
the origin of possible ferromagnetic instability for a (001) 
LAO-STO interface and then apply the results to other 
interfaces and materials. The main result of this paper 
is that, despite the aforementioned results for the single- 
band Hubbard model, we find that the directional char- 
acter of the t2 g orbitals allows for the possibility of weak 
ferromagnetism in this multi-band (or multi-orbital) sys- 
tem. Such ferromagnetism may occur precisely in this 
intermediate-high carrier density regime, and is strongly 
dependent upon the microscopic details such as orbital 
splittings and the magnitude of inter-orbital Hund's cou- 
pling. 

For our discussion, we will require a few particulars 
of the conduction band states in STO, which are well- 
established. The low lying octahedral t2 g crystal field 
levels of Ti comprise yz, xz, and xy orbitals. Owing 
to its directionality, hopping amongst xy orbitals occurs 
most strongly with neighbors in the x-y plane, and much 
more weakly in the z direction (a relative strength of 
approximately 1/8 is assumed for the ratio of z axis to 
x-y hopping in Ref. H2I ) ; the same considerations apply 
to yz and xz orbitals with appropriate rotations. Thus 
in a bulk system with cubic symmetry, there are three 
bands each of which disperses predominantly in 2 of the 
3 cartesian directions. When confinement is introduced 
in the z direction for a (001) inteface, the states become 
quantized and bound in this direction, leading to two 
dimensional subbands. The xz and yz bands, which dis- 
perse in the z direction, suffer the greatest loss in kinetic 
energy, and are lifted somewhat in energy. The lowest 
energy subband is therefore expected to be of xy char- 
acter, while in higher subbands the xz and yz states ap- 
pear. The lowest xy subband disperses fairly uniformly 
in the two dimensional plane, while the higher xz and yz 
subbands become approximately one- dimensional, owing 
to the quenching of one of two of their main dispersion 
directions. 



The reduction of the kinetic energy of the xz and 
yz subbands suggests we consider them for possible fer- 
romagnetic polarization. We therefore adopt a mini- 
mal model with three two-dimensional sub-bands for the 
(001) interface, with the Hamiltonian H = TLq + Hi, 
where the kinetic energy is 
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Here do a ,d xa ,d ya describe the xy, xz, and yz bands 
(with spin polarization a —^,],), respectively, m, is an 
effective mass, and A is the subband crystal field split- 
ting. We have assumed a tetragonal crystal field symme- 
try for the (001) interface in Eq. ([1]), so m x = m y . Since 
we always consider the band bottom, this is equivalent 
to taking a tight binding model with hopping amplitude 
ti = l/(2m,;a 2 ), where a is the lattice spacing. Here we 
make an approximation to « t x = t y = t so that all the 
effective masses are equal. We take on-site interactions, 
of the form 

Hi = J7^n lt (r)n. 4 (r) + U' £ ^(rjn^r) 

r.i r ^7^j 

-J H Y i S i {T)-S i {v), (2) 

where n ia (v) = d\ a (r)d ia (r), n»(r) = Y, a n "x( r )> and 
s *( r ) = lEa^lW^^W- As usual, we expect 
the intra-orbital interaction Hubbard U to be the largest 
interaction, with the inter-orbital interaction U' and the 
Hund's coupling Jh rather smaller, U'/U, Jh/U < 0.3. 

To analyze this Hamiltonian, we treat the U' and Jh 
as small (which they are, relative to U), and consider 
possible ferromagnetic instabilities they induce. With 
U' = Jh = 0, the Hamiltonian is decoupled to three 
single-band problems, and hence does not support fer- 
romagnetism. What we require for the analysis is the 
susceptibility of each orbital subsystem, including the ef- 
fects of strong on-site U. For the xy subband, which 
is two-dimensional, an on-site interaction U has little ef- 
fect, and in the low density (per lattice site) limit studied 
here the interactions can be exactly treated by a stan- 
dard T-matrix ladder summation. The result is simply 
a Fermi liquid with small Landau parameters, which can 
be neglected at the level of the present consideration. 
Neglecting the Fermi liquid correction, the spin suscepti- 
bility of the xy band per site is then X2d = ma/ (27r) = 
l/(47r£). For the xz and yz subbands, however, being 
one-dimensional, the effect of U is critical. As is well- 
known, an arbitrary strength repulsive interaction in one 
dimension transforms a Fermi liquid into a Luttinger liq- 
uid. The Luttinger liquid is described at low energy 
by bosonization, which decomposes the Hamiltonian of 
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the one dimensional system into charge and spin modes, 
with independent propagation velocities v c and v s , re- 
spectively, which generalize the Fermi velocity vp. No- 
tably, when the fractional electron occupation x\ of a 
one-dimensional chain is small, < 1, the effect of in- 
teractions is strongly enhanced (this is a result of the 
fact that, in one dimension, an arbitrary repulsive scat- 
terer becomes perfectly reflecting). As a consequence, 
the spin velocity is strongly suppressed. In the large U/t 
limit this can be obtained from the results of Shiba|19j for 
the susceptibility xid, and the Luttinger liquid relation 
Xid = a/2irv s , which gives 



Jax 1 , 



(3) 



where J = ir 3 t 2 /(3f7), modeled as a chain with hopping 
t. We expect that Eq. ([3]) holds for arbitrary U/t and 
x\ -C 1, but with J becoming of order U for U -C t. 
Notably, in the same limit, the charge velocity is much 
larger, v c ~ tax\. The smallness of the spin velocity 
indicates the strong suppression of spin exchange in this 
limit, and as a consequence a much reduced energy cost 
to polarize the electrons in the xz and yz bands. 

Consider now the energy change of the system on po- 
larizing the electrons. If the spin (per site) in the xy 
band is Mq and that in the xz and yz bands is Mi, this 
is (per site) 
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M l + o M i - Jh(2M M 1 + Ml) (4) 

The ground state is paramagnetic if this quadratic form 
is positive semi-definite, as it is for a small Jh- Con- 
versely, a ferromagnetic solution with Mq , Mi ^ is fa- 
vored when it becomes unbounded below, which occurs 
when Jff(Jff + l/(2x2d)) > l/( 2 XidX2d)- Using the forms 
of the one-dimensional and two-dimensional susceptibili- 
ties, this condition gives 
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This instability means Mo, Mi will increase until Eq. Q 
becomes invalid and stabilize the system, and the natural 
guess is that this occurs when the less occupied excited 
subbands become fully polarized. Remarkably, this oc- 
curs for arbitrarily weak Hund's coupling Jh, provided 
the filling of the upper xz and yz subbands is sufficiently 
small (see Eq. (^J), and of course non-zero. This gives 
a mechanism for fcrromagnctism at intermediate carrier 
density, when the total density is near the critical value 
needed to just populate the xz and yz subbands, with 
magnetism disappearing both for smaller and larger car- 
rier density. 

It may appear that Luttinger liquid physics is essen- 
tial to this above picture, a scenario that might seem 
overly exotic and restrictive. However, this is not the 
case. What is really required is simply an enhancement 



of the magnetic susceptibility of the upper xz and yz 
subbands when their occupation is small. The same ten- 
dency should be present even in the presence of weak 
dispersion normal to the direction of orbital extent, e.g. 
in the y direction for the xz orbital. This effect may 
be included in the model by adding a small t' hopping 
term, and would be expected to cause a crossover from 
Luttinger liquid behavior to Fermi liquid behavior at low 
energy. By continuity, for small t' jt, this Fermi liquid 
must have an enhanced spin susceptibility captured by a 
large Fermi liquid correction Fq. However, the eventual 
two-dimensionality induced by non-zero t' must prevent 
the divergence of the spin susceptibility of the xz and 
yz bands at small Xi, and if this effect is too large, the 
ferromagnetic instability may be entirely removed. To es- 
timate this, we note that the one-dimensional magnetic 
susceptibility divergence is cut off when the Fermi en- 
ergy of the xz and yz bands (measured from their min- 
imum) is comparable to the hopping t' , i.e. tx\ ~ t' . 
The same condition describes the change from an open 
Fermi surface to an elliptical one. This gives the condi- 
tion if < Jut I J for the ferromagnetic phase to occur (we 
neglect numerical prefactors here due to the imprecision 
of the matching argument). 



Interfaces 


Symmetry 


Local orbitals 


(001) 


4- fold rotation 


xz, yz; xy 


(110) 


2- fold rotation 


-^(xz + yz); xy; ^(k - yz) 






-js(xy + e l ~%-yz + e~ l2 S~xz), 


(111) 


3- fold rotation 


-^(xy + e~ i2 fyz + xz); 
75 (xy + yz + xz) 



TABLE I. The relevant local crystal symmetries and local 
orbital states for different interfaces. ';' delimits the sets of 
locally degenerate orbital states. Note that there could be a 
small hybridization between (xz + yz) and xy orbitals for 
the (110) interface. 

For the (001) interface, it is the tetragonal crystal 
field environment and the quasi-two-dimensional geom- 
etry of TiC>2 layers that lead to the particular orbital 
and band structure discussed above. Recently the (110) 
and (111) LAO-STO interfaces have also been prepared 
experimentally 2^, 21 1 , and both interfaces appear to sup- 
port STO electron gases, though the (111) interface is 
polar and the (110) is not. As listed in Tab.|H these two 
interfaces have different local crystal field environments 
and hence different local orbital configuration from the 
(001) interface. Can these two interfaces also support 
fcrromagnctism - under ideal disorder-free conditions - 
at certain electron fillings? 

As usual, the e 9 doublets are always higher in en- 
ergy and do not play any role. For the (110) inter- 
face, the three t2 S orbitals are splitted into three non- 
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degenerate orbitals, -^( xz + yz), -^(iz — yz) and xy. 
In the first approximation, the local hybridization be- 
tween -^(xz + yz) and xy orbitals may be neglected. 
When these two orbitals form bands, they are also quasi- 
one-dimensional just like xz and yz orbitals for the (001) 
interface. Hopping among -^(xz + yz) (/xy) orbitals oc- 
curs most strongly with neighbors along z ( /[110]) lattice 
directions. The ^(xz — yz) subband is two dimensional 
and its band bottom is the lowest among the three sub- 
bands. Due to the reduced symmetry of the (110) inter- 
face, the two upper quasi-one-dimensional subbands are 
split. Based on our above discussion of ferromagnetic in- 
stability for the (001) interface, we also expect emergent 
ferromagnctism for the (110) interface when the filling of 
the quasi-one-dimcnsional subband is sufficiently small. 
Because the two upper quasi-one-dimensional subbands 
are not degenerate, there may even exist two ferromag- 
netic regime as the electron filling of the two upper sub- 
band is increased. One should note that the discussion 
here assumes no hybridization between -^(xz + y z ) an d 
xy orbitals. In reality, there are always small hybridiza- 
tion between these two orbitals. If this hybridization is 
very small (smaller than 0(y/ Jh / J)), the resulting two- 
dimensional Fermi liquid should still have a large spin 
susceptibility and ferromagnctism can still be present. 

For the (111) interface, although locally the crys- 
tal field splits three t 2g orbitals into one &i g state, 
^(xy + yz + xz), and two e' 2g states, -^{xy + e l ~^~yz + 
e~ z ^~xz) and -^(xy + e~ l ^~yz + e l ^~xz), the electron 
hopping strongly hybridizes three orbitals and leads to 
two-dimensional Fermi liquids. Therefore, our mecha- 
nism for ferromagnctism cannot predict a ferromagnetic 
state. 

In the above, we have discussed possible itinerant ferro- 
magnctism in two dimensional STO systems. Other the- 
oretical works have instead proposed mechanisms relying 
on localized electron moments. While we believe that 
Mott localization of electrons near the interface should 
not occur for ideal structures, sufficient disorder and in- 
teractions together might create some truly localized mo- 
ments. If the localized electron mechanisms are correct, 
we predict significant dependence of the ferromagnctism 
on disorder, and indeed that it should weaken as sample 
quality is improved. The itinerant mechanism discussed 
here has its own distinct predictions. For example, vary- 
ing the electron concentration away from the critical den- 
sity by tuning the back gate voltage may easily suppress 
the ferromagnetism, and ferromagnetism should be ab- 
sent at the (111) interface; neither prediction applies for 
the local moment mechanism j22j| . 
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